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A pulse vaccination SIR model with periodic infection rate (i{t) have been proposed and studied. 
The basic reproductive number Rq is defined. The dynamical behaviors of the model are analyzed 
with the help of persistence, bifurcation and global stability. It has been shown that the infection- 
free periodic solution is globally stable provided Ro < 1 and is unstable if i?o > 1. Standard 
. . . bifurcation theory have been used to show the existence of the positive periodic solution for the case 

of -Ro — > 1^. Finally, the numerical simulations have been performed to show the uniqueness and 
, the global stability of the positive periodic solution of the system. 
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(-H ' Transmissible diseases have tremendous influence in human life. Every year billions of people suffer or die in 
various infectious diseases. In recent times, the emerging and reemerging communicable diseases have led to a revive 
interest in the study of infectious disease. Mathematical models are widely used to understand the mechanisms of 
spread of infectious diseases, and become an important tool to analyze the spreading and controlling the diseases 
[1]. Epidemiology modelling and mathematical study on those model can contribute to the design and analysis of 
• epidemiological surveys, suggest some crucial conclusion that should be collected. By determining the key parameters, 
, and finding it's effects of changes in the parameter values, it identify the trends, make general forecasts, and finally 
' estimate the uncertainty in forecasts [2]. Another important motivation of the development of this field is the 
I evaluation of various vaccination/control strategies for human as well as animal. Vaccination is an effective way 
to control the transmission of a disease among the species. Mathematical modelling can contribute to the design 
and assessment of the vaccination strategies. A theoretical examination of pulse vaccination policy for SIR model 
. have been shown by Stone et. al [3]. They found a disease- free periodic solution with the same period as the pulse 
vaccination, and studied the local stability of this solution. Nokes and Swinton have studied the control of childhood 
viral infections by pulse vaccination [4]. Major application of the pulse vaccination method for SIR and SEIR epidemic 
models have been studied by A.d'onofrio and his group [5-7,9]. They have studied global asymptotic eradication in 
the presence of vaccine failure [5]. Yicang zhou and Hanwu Liu studied the stability of Periodic Solutions for an 
■ SIS Model with Pulse Vaccination [8]. Fuhrman et. al studied asymptotic behavior of an SI Epidemic Model with 
Pulse Removal [10]. The global stability of the disease-free periodic solution for SIR and SIRS models with Pulse 
Vaccination have been investigated by Jin [11]. 
' Recently, pulse vaccination, the repeated application of vaccine over a defined age range, is gained prominence as 
a strategy for the elimination of childhood viral infectious such as measles and poliomyelitis. This policy is based on 
the suggestion that measles epidemics can be more efficiently controlled when the natural temporal process of the 
epidemics is antagonized by another temporal process [12,13]. 

On the other hand , in the last 30 years epidemiology of infectious diseases has focused on the theoretical and the 
numerical-experimental study of the effects of periodically varying external factors on the time course of the incidence 
of infectious diseases. It is well known that the epidemic models with constant contact rate are not able to reproduce 
realistic endemic situations, since incidence of many infectious diseases are not constant [14]. The endemic nature of 
many communicable diseases is characterized by a wide range of temporal oscillatory patterns: annual or polyannual 
periodicity [15], apparently random oscillations with a strong annual component in the power spectrum [16] etc. 
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In this paper, we have studied the dynamical behavior of a pulse vaccination SIR model with periodic infection rate. 
The paper is organized as follows: Section 2 gives the SIR model with pulse vaccination and definition of the basic 
reproductive number i?o, and demonstrates the existence of an infection-free periodic solution. The global stability of 
the infection-free periodic solution is obtained in Section 3. In Section 4, we discuss the persistence of the infectious 
disease by obtaining the uniform persistence conditions i?o > 1 of the infectious disease. Section 5, concentrates 
on the existence of the positive periodic solution of the periodic SIR model with pulse vaccination. The existence 
of positive periodic solution is obtained by using a well-known bifurcation result of Rabinowitz [17] for Rq — > 1+. 
Finally, numerical simulations and a brief discussion conclude the paper. 



Our SIR model is based on the following assumptions: 

1. The total population size is constant and is denoted by N, with N ^ S + I + R, the population is divided into 
three groups: (a) The susceptible class, S, comprising those people who are capable of catching the disease; (b) The 
infectives, I, comprising those who are infected and capable of transmitting the disease; (c) The recovered class, R, 
comprising those individuals who are immune. 

2. The total population size is N and the per capita birth rate is a constant b. As births balance deaths we must 
have that the per capita death rate is also b. 

3. The population is uniform and mixes homogeneously. 

4. The infection rate (3{t) is defined as the total rate at which potentially infectious contacts occur between two 
individuals. A potentially infectious contact is one which will transmit the disease if one individual is susceptible and 
the other is infectious, so the total rate at which susceptibles become infected is (3{t)SI: Biological considerations 
mean that (3(t) is a continuous function. We also assume that f3(t) is not identically zero, positive, non-constant and 
periodic function with period uj > 0. 

5. The infectives move from the infectives class to the recovered class at a constant rate 7 where (I/7 ) is the 
average infectious period conditional on survival to the end of it. 

6. We also assume that a is a non-negative constant and represents the death rate due to the disease (that is, the 
disease-related death rate). 

On the basis of the above assumptions a SIR model can be written as a set of coupled non-linear ordinary differential 
equations as follows: 



Here the parameters K, b, 7 are all positive constants. The constant bK = A is the inflow rates, b is the per capita 
birth rate and the per capita death rate, so that K represents a carrying capacity, or maximum possible population 
size. 

Adding the three equation of system (2.1), we get 



For small value of the parameter a the total population size become closer to the carrying capacity K. 

According to usual convention, pulse vaccination can be defined as the repeated application of vaccine across an 
age range. Let us assume the pulse scheme proposes to vaccinate a fraction of the entire susceptible population in 
a single pulse, applied k time in time period [0, w]. Pulse vaccination gives lifelong immunity to p^-S susceptibles who 
are, as a consequence, transferred to the recovered class (i?) of the population. 

When pulse vaccination is incorporated to SIR model (2.1), the model takes the form as follows: 



2. The SIR model without and With pulse vaccination 



S' ^bK~ (3{t)SI - bS, 
r =l3{t)SI ~{b + -f + a)I, 
R' = jl - bR. 



(2.1) 



N'{t) =bK -bN - al. 



(2.2) 




(2.3) 




(2.4) 
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where < Pk < G Z+) are constants and g > is an integer such that Pk+q = Pk,Tk+q = Tk + oj. Note that the 
dynamics of the total population size N{t) still satisfy equation (2.2). 
From biological view point, we can assume that the domain 

{{S,I,R) : S >0,I >0,R>0,S + I + R < K} 

is a positive invariant set of system (2.3) and (2.4). 

Let us start to analyze system (2.3) and (2.4) by demonstrating the existence of an infection-free periodic solution, 
in which infectious individuals are entirely absent from the population permanently i.e., 

/(t)=0,t>0. (2.5) 

Under this conditions, the systems (2.3) and (2.4) reduce to 

' S' = bK-bS, t^Tk. 
R'^-bR, fc = 0,l,2-... 

S{T+) = {l-pk)S{rk), ^'■^> 

[ R{T+)=R{Tk)+PkS{Tk). 

and N' = bK — bN , therefore N{t) K &s t ^ oo. For convenient, we take N = S + R = K. Then, solving the 
equation (2.6), we get 

S(t)=W{t,0)SiO) + bK [ W{t,T)dT 

Jo 

and 



R{t) = W{t, Q)R{Q) + ^o<r,<tW{t, Tk)pk, 



where 



W{t,r)= n (1-P,)e- 



,-6(t-r) 

T<Tj<t 



Since W{uj,0) — Yii^ ^ Pj)^ < 1; therefore equation (2.6) has a unique w-periodic solution (S'*(<), 0, i?*(t)) with 
the initial conditions S*{0) = bK W{uj,T)dT/{l - 1^(^,0)), i?*(0) = J2 l^(w,rfe)/(l - W{uj,0)). 

0<Tk<t 

Let us define the basic reproductive rate of model (2.3) and (2.4) as follows: 

uj{b + Q! + 7) 

where S*(t) is the periodic infection-free solution. 

3. Local and global asymptotic stability of the infection free solution 

In this section, we will prove the local and global asymptotic stability of the infection free solution {S* (t), 0, R*{t)). 

The local stability of the w-period solution {S* {t),0, R* (t)) may be determined by considering the linearized SIR 
equation of (2.3), (2.4) about the known periodic solution (5*(t), 0, i?*(t)) by setting S{t) = S*{t)+x{t), I{t) = y{t), 
R(t) = R*{t) + z(t), where x{t),y{t), and z(t) are small perturbation. The linearized equations may be written as 

x{0) 

m I y{0) 

m 




where $(t) = (pij{t),i,j = 1,2,3 satisfies 

d$(t) _ 
dt ^ 



-b -[3{t)S*{t) 

(3{t)S*{t) - (6 + 7 + a) I $(t), 
7 -6 
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with $(0) = E, where E is the identity matrix. Therefore, the system ( 2.4) become 

I-Pk 0\ f x{Tk) 
10 y{n) 

Pk I J \ z(Tfe) 

Hence, according to the Floquet theory, if aU eigenvalues of 

/ "J \ 

/ Ua-Pk) \ 

k=l 

1 

. n(i-Pfc) 1 

\ fc=2 / 



M 



$(tj) 



are less than one, then the w-periodic solution {S*{t), 0, R*{t)) is locally stable. After calculation we get. 



where 



"^322 
^32 



Vl2 

^{t) = I ^22 

^32 e-" 



exp(/o*[/3(r)5*(r)-(6 + a + 7]dT) 

- exp i-bt) /3(T)5*(T)(^22(T)e^-dT, 

exp —bt 7(^22 {T)e'''^dT. 



Thus, if ^1, ^2 and /is are the eigenvalues of matrix M, then they are given by, 

q 



Ml -Pk)e '"^,fJ,3 



k=l 



and 



/i2 = f22i^^) = cxp( / (3{T)S*{T)dT - (6 + a + 7)ti;. 
Jo 



Obviously /^i < 1, /is < 1 and /t2 < 1 if and only if 

/;/3(r)5*(r)dT 



a; 



< 5 + a + 7. 



Therefore, the periodic infection-free solution (S'*(i), 0, i?*(i)) is asymptotically stable if i?o < 1. Now we are in a 
position to summarize the above results in the following theorem: 

Theorem 3.1. The periodic infection-free solution {S*{t), 0, R*{t)) of system (2.3) and (2.4) is local asymptotically 
stable provided i?o < 1. 

In order to prove the global stability of infection-free solution {S*{t), 0, R*{t)), we need the following lemma. 
Lemma 3.1. (Comparison theory [18]). Assume that m £ PC[R+,R] with points of discontinuity at < = t„ and 
is left continuous at t = t„, n = 1, 2, • • • , and 



D_m(t) < g{t, m(t)), i ^ r„, n = 1, 2, • • ■ , 
"^(t+) < ll)n{m{Tn)),t = T„,rt = 1,2, ••• , 

where g S C[R+ x ■)/)„ S C[R,R\ and tpniu) is nondecreasing in u for each n = 1,2, 

maximal solution of the scalar impulsive differential equation 



(3.1) 

Let r{t) be the 



u' = g{t,m{t)),t ^ T„,n^ 1,2, ■• ■ , 
m(t+) = Tpn{u{Tn)),t = r„ , n = 1, 2, • • 

M(t+) = Mo, 



(3.2) 
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existing on [to,oo). Then < uq impHes m{t) < r{t). 

Remark. In Lemma 3.1, assume the inequahties (3.1) reversed. Let p{t) be the minimal solution of (3.2) existing 
on [to, oo) . Then, m{t) > p{t). 

Lemma 3.2. For the following impulsive equation 

S' = bK-bS, t^Tk. „^ 



SiT^)^{l-Pk)S{n), n = l,2,... , 

has a unique positive uj -periodic solution S*{t, Sq) for which 5*(0, 5*0) = 5*0 and S*{t, Sq) > 0, t G R+, and S*{t, Sq) 
is global asymptotically stable in the sense that lim \S{t, So) — S*{t, Sq)\ = 0, where S{t, So) is any solution of system 

t — >oo 

(3.3) with positive initial value S{0, Sq) = Sq > 0. 

Proof. The first part of the proof follows from the results of Section 2. Therefore we have to show only 
lim \Sit,So)- S*{t,S^)\ =0. 

t — >-oo 

Since 

\S{t,So)- S*{t,S*Q)\^W{t,0)\SQ- S*o\, 
thus we need to show W{t, 0) — > as ^ 00. Suppose t e {ntu, {n + l)uj], then 



wit,o) = n (i-p.) 



-6t 



^ ^ e 

0<Tj<t 

= n (1 -Pi)e~''"'^ n (1 -Pj)e-''(*-"") 

< Tj < n nuJ<Tj<t 

= [ n (l ~ Pj)e~'"^Y' n (1 -Pi)e-''(*-"'^' 

0<Tj <.UJ 7lU)<rj <t 

< n n (i-p.K 

0<Tj<uj 0<rj<uj 



-boj]n 



Thus hm W{t, 0) = 0, since J] (1 ~ P])^^ < 1- Hence the proof. 

Now, we are totally ready to show the global asymptotic stability conditions of the infection free solution. We claim 
the following theorem: 

Theorem 3.2. The periodic infection-free solution {S*{t),0, R*{t)) of system (2.3) and (2.4) is globally stable if 
Ro < 1. 

Proof. From the first equation of (2.3) and (2.4), and by using Lemma 3.1 and 3.2, we obtain that, for any given 
< e < u}{l - RQ){b + 'f + a)/2 (3{t)dt, there exists Ti > 0, such that 

S{t) < S*{t)+e, Vt>Ti. (3.4) 

By substituting (3.4) into second equation of (2.3), we obtain 

= Pit)SI ^{b + -f + a)I< P{t)S*{t)I{t) - (6 + 7 + a)I{t) + eP{t)I{t) 

dt 

By comparison theory, for t G (Ti + noj, Ti + (n + 1)0-"], we have 

I{t) < I{T^) exp(/* [/5(r)5* (r) - (5 + 7 + a) + P{T)e])dT 

= liTi) exp{/|^+""[/3(T)5*(r) - (6 + 7 + a)]dT + /3(r)dr 

+s Ik+n. Pi^)dr + /t,+„J/3(^)5* (r) - (6 + 7 + «)]rfr} (3.5) 

< /(Ti) exp{n !^[I3{t)S*(t) - (5 + 7 + a)]dT + (n + l)e /q" /3(r)dT + /3(T)5*(T)dT 

< B exp{nw(6 + -f + a){Ro - I) + ne ;fl(T)dr} 

< B cxp{nuj{b + J + a)^{Ro - 1)} 

where B = /(Ti) cxp{ /?(r)(£ + S*{t))dT}. Thus if i?o < 1, it follows from the above equation that lim I{t) = 0. 
Let H{t) = \S{t) - S*{t)\, then we have 

D+H{t) = sign(5(t) - S*{t)){S'{t) - S*'{t)) < -bH{t) + (3{t)SI, (3.6) 
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H{T+) = {l-pk)H{n). (3.7) 

Since S(t) < K, f3(t) < max /3(t) = f3* , the equation (3.6) reduces to 
te[o,u] 

D+H{t) < -bH{t) + KI3*I{t). (3.8) 

For t G {Ti + TiLo, Ti + (n + l)cj] , from (3.5), we see that I[t) < Se""", where 9 = ^Lu{b + j + a){l - Rq). By lemma 
3.1 and (3.7), (3.8), we have 

H{t)<H{T^) W (l-p,)e"''(*-^i r J| (1 - j9,)e-^(*-^)/(r)dT. (3.9) 

Ti<Tj<t "^^1 T<Tj<t 

Since 



= a - P,)e-^(*-)/(r)rfr + fj^^^^^ ^Jljl p,)e-'i'-^) I{r)dT + ■■■ 

+ iirtSTr a - P,)e-^(*-)/(r)dr + fU^^^Jl - p,>-''(*-)/(r)dr 

< ^UJl Pi)e-^<~'-^)dT + ^ a - P,)e-''(*-)e-«dr + . . . (3.10) 
+ /(:-tS+Trn (l-p.)e-''(*-)e-("-i)«dr+}£^^ H (1 - P,)e-^(*-)e-'^"dr] 

^ ' T<Tj<t T<Tj<t 

= Be^^i(e'"^ - l)e-''*[r?" - e"^'"^-^)]/^!'? ^ ef*"^-^)) 

< B(e*" - l)e-''""[77" - e"(''"-^)]/6(?/ - e^*""-*)), 

where < r; H (1 -Pj) < 1- From (3.9) and (3.10), we obtain that hm H{t) = 0, i.e., S{t) S*{t) as i ^ cx). 

j=l t^oo 

Similarly, we can prove R{t) R*{t) as i ^ cx). The proof is completed. 

4. The uniform persistence of the infectious disease 

In this section, we will discuss the uniform persistence of the infectious disease, that is, lim inf I{t) > a if i?o > 1. 

t — >QO 

To discus the uniform persistence, we need a lemma first. 

Lemma 4.1. If Rq > 1, then the disease uniformly weakly persists in the population, in the sense that we will 
able to find a constant c > such that lim sup/(i) > c for all solutions of (2.3) and (2.4). 

t — >oo 

Proof. Let us suppose that for a given e > 0, there exists a solution with lim sup/(i) < e. From the first equation 

t — >oo 

of (2.3), we have 

S' = bK - l3{t)SI -bS> bK - p*Ke -bS,t^Tk. 
Consider the following equation 

u' ^bK - P*Ke -bu,t^Tk. , . 

u{T+)^{l-pk)u{T,),k^l,2,--- ^^■'> 

By lemma 3.2., the (4.1) has a unique positive uj -periodic solution u*{t) for which u*(0) = Uq and u*{t) is global 
asymptotically stable. So 



Let 



W{t,0) [!:'W(uj,T)dT /•* , 
A = f3*K max { ^ ' ' ^ + / W{t, r dr}. 
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By (4.2), we see that 

u*{t)>S*{t)-Ae. (4.3) 

By comparison theory, we obtain that 

>m[mu{t)-{b+^+a)]. ^^■^> 

Since, u*{t) is global asymptotically stable, for our previous e above, there exist Ti > 0, such that u{t) > u*{t) — e, 
t > Ti. From (4.3) and (4.4), we get 

l'it)>I{t)[(3{t)S*it)^{b + -f + a)-eil + A)P{t]. (4.5) 

Integrating over intervals [Ti , t] we obtain 

I{t) > exp(/* [/3(t)5*(t) - (6 + 7 + a) - e(l + A)/3(t)]dT 

= /(Ti) exp{4^+"n/3(r)^* (r) - (5 + 7 + a)]dr - £(1 + A)l3{t)nu ,^ 
+ Sk+nM^)S*{r) - (6 + 7 + a)]dr - e(l + A)/3(t)(t - nc.)} 
> Cexp{nu;(6 + 7 + a)(i?o - 1) - e(l + A)/3*nw}. 

where t £ (Ti + nw, Ti + (?i + l)uj], C = /(Ti) exp[-(6 + 7 + a + ^(1 + A)/3*)w]. Taking 

(6 + 7 + a)(i?o-l) 



< e < 



2/3*(l + A) 



thus I{t) — > 00 as t — > 00, a contradiction to the fact that I{t) is bounded. This finishes the proof. 

Theorem 4.1. If i?o > 1, then the disease uniformly persistence, it is that there exists a positive constant tr such 
that for every positive solution of (2.3) and (2.4), 

lim inf /(t) > a > 0. 

t — ^00 

Proof. Let 

bK{R^ - 1) 



< ?7 < 



2i?o(a + 7)' 



It can be obtained from the Lemma 4.1 that for any positive solution of (2.3) and (2.4) there exists at least one to > 
such that /(to) > ?7 > 0. Then, we are left to consider two possibilities. The first case is I{t) > ?/ for all large t > to. 
The second one is I{t) oscillates about 77 for large t. The conclusion of Theorem 4.1 is obvious in the first case since 
we can choose a = rj. 

For the second case, let ti > to and t2 > ti satisfy 

/(ii) = /(fj) = ?7, and I{t)<i] for ti<t<t2. 

Next, we introduce a new variable V = S + I, and it follows from the first two equations of (2.3) and (2.4) that 

V = hK^bV-{^ + a)I, t^Tk. 

/' =/3(t)(V^-/)/- (6 + 7 + a)/, fc = 0,l,2---. ^^■'> 



If I{t) < T], then 



ViT+) = {l~Pk)ViTk)+PkIiTk), 
I{T+)=I{Tk), fc = 0,l,2. 



V = hK-bV -{^ + a)I > (1 + _ t^ Tu. 
y(r+)>(l-Pfc)y(r,)). 



(4.8) 



(4.9) 
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Consider the following equation 



X' 



I _ bK 



By using Lemma 3.2, we can obtain equation (4.10) has a unique positive uj- periodic solution x*{t), and x*{t) is 
global asymptotically stable in the sense that lim \x{t, xq) —x*{t)\ = 0, where x{t, xq) is any solution of system (4.10) 

t — ^oo 

with positive initial value x(0, xq) = xq > 0. By (3.3) and (4.13) we can easily get 

A Ho 

where S*{t) is a unique positive u- periodic solution of system (3.3). 

Then the comparison principle and the global asymptotically stable of x*(t) implies that there exists a positive 
constant Ti > 0, such that 

Vit) > ^{1 + ^)S*{t), for all t>ti+Ti. (4.11) 
For above ti and rj the solution of the initial value problem 

y'{t)^P{t)[Ul + ^)S*{t)-y]y-ib + -/ + a)y, y{h) = 7^ (4.12) 

is 

^ V^MI', m + ^mr)S*{r) - (6 + 7 + a)]dr) 

^ 1 + ,,/* /?(r)exp(/,;[i(l + ^)/3(0)5*(0)-(6 + 7 + «)Wdr 

The differential equation in (4.12) has a periodic solution 

^ y* ^Mll m + i;mr)S*ir) - (& + 7 + a)]dr) 

^ 1 + r /3(r) exp(/,;[i(l + i^)mS*{e) - (6 + 7 + a)]d0)dT 

where y = —rmz ; — , — ; • omce the condition /tn > 1 of ihcorcm 4.1 implies that 



/o"[i(l + -^)/3(T)5*(r) - (6 + 7 + a)]dT = -^(b+j+aKB.o-i) ^ ^^^^ ^-^-^^ 



lim exp( / [i(l + -L)/3(r)5*(T) - (6 + 7 + a)]dT) = oo, 

lim / /3(T)exp(/ [i(l + -^)/3(0)5*(0)-(6 + 7 + a)]rf0)rfT = oo 
Jti 2 

and the expression of y(i) — y*{t) it follows that 

lim |y(i)-y*(OI =0. (4.15) 

r — >oo 

From (4.15) we see that there exists a positive constant T2 > such that 

y{t) > p=\ min y*(t) > 0, for ah t > ti + T2. (4.I6) 

Let T = maxjTi, T2} and define 

(J = min{p, ryexp(— (6 + 7 + a)T}. 
If t2 — ti < T, from the following equation 



I' = p{t){V - /)/ - (6 + 7 + a)I, 



(4.17) 
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we have the inequaHty 

I'{t) > -{b + -f + a)I, 



and the comparison principle implies that I{t) > T]exp{ — {b + 7 + a)} > T^expl — (fe + 7 + a)T}, i.e., I{t) > a for all 
t G {t\,t'2). If t2 — ti > T, we divide the interval [^1,^2] into two subintcrvals [ti,ti +T] and [ti +T,t2]- I{t) > <t is 
obvious in the interval [ti, ti +T]. By (4.11) and (4.17) we see that in the interval [ti + T, ^2], 

I'it) > P{t)[hl + ^)S*it) - /]/ - (5 + 7 + a)I. (4.18) 

Therefore again by using comparison principle we get I{t) > y{t) > p > cr for t E [ti -\- T, t2]- The above analysis 
is independent of any interval [ii,t2]! and the choice of a is independent of any positive solution of (2.3) and (2.4). 
Therefore, the persistence is uniform for all positive solution. 



5. Existence of the positive a;-periodic solution and bifurcation 



Let PC{J,R) {J C R) be the set of continuous functions ip : J ^ R ior t. E J,t ^ Tk, that have discontinuities of 
the first kind at the points Tk E J where they are continuous from the left. Let PC'{J,R) be the set of functions 
ip : J R with derivative ^ G PC{J,R). We consider the Banach spaces of w-periodic functions PCuj = {"0 G 
PC([0,a;], i?)| V(0) = 0(w)} under the supremum norm || -0 ||pc^= sup{| if) \: t € [0,a;]}, and PC'^ = {-0 S 
PC"([0, w], i?)[ ■0(0) = 0(w)} under the supremum norm |[ tj: |Ipc^= max{|| tp \\pCuji II 0' IIpc^}- We will also consider 
the product space PCui x PGu which is also a Banach space under the norm || (0i, ■0i) ||pc^ = || ^/'i ||pc^ + || 02 ||pc^- 
Moreover, for any f & C^j (or PC^j) we define the average of / by / := f{s)ds/uj. 

We need the following lemmas which are about linear impulsive periodic equations (also see [20,21] ). 

Lemma 5.1. Suppose Cij{t) S PC^j, and 
q 

(1). if cii 7^ ;^ ln[ n rr^] ^^'^ ^22 7^ 0, then the linear homogeneous periodic impulsive equation 



k=l 



X[ = CiiXi{t) + Ci2X2it), t^Tk 
x'2 = C22X2{t) , /s = 0,l,2---. 

^li^k) = - Pk)xi{Tk), t = Tk 

has no nontrivial solution in PC^ x PCoj. In this case, the nonhomogeneous system 

y'l = Cliyi{t) + Ci2?/2(t) + fl, Tk 
y2 = C222/2(t)+./2, fc = 0,l,2---. 

yii'^k) = i'^ - Pk)yiin), t = Tk 
y2{T^) = y2{Tk), 



(5.1) 



(5.2) 



has a unique solution (2/1,2/2) £ PC^^ x PC^ for every (/i,/2) S PC^ x PC^^ and the operator L : PC^ x PC^^ 

PC^ X PC^ defined by (j/i, 1/2) = ^(/i, /2) is linear and relatively compact. 
q 

(2). if cii 7^ — ln[ -T-^ — ] and C22 — 0, then (5.1) has exactly one independent solution in PC^j x PC^. 

" fc=i 
Proof. (1) Since 

2^2 (0 = 2^2 (0) exp / C22{s)ds (5.3) 
Jo 

the condition C22 7^ implies that X2 ^ PC^i unless X2 = 0. Then xi{t) = Yl — Pfc)a^i (0) exp cii{s)ds and 

0<Tk<t 

q 

cii ^ — ln[ n rr — ] in turn implies x'l ^ PCui unless xi = 0. 
" k=i 
In this case, 



2/2 = C22y2{t) + f2 
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has a unique solution j/2(i) G PC^ and the operator L2 : PC^i PC^ defined by 1/2 = ^2/ is Hnear and relatively 
compact. Furthermore, 



viijk) = {'^ - Pk)yi{Tk), t^Tk, 



foi' /s — C12-ZJ2/2 + /i G ^'C'w has a unique solution (since cn ^ ^ H ]) hi PC^j and yi = L1/3 defines a linear, 

fc=i ' 

relatively compact operator 

Li : PC'i^ > PCi^. 

Thus (5.2) has a unique w-period solution in PCi^ x PC^^ given by {yi, j/2) = L{fi, j'l). where 

i(/l,/2) - (il(ci2i2/2) + /l),i2/2). (5.4) 

(2). Under the stated assumptions, X2 as given in (5.3) lies in PC^ for all initial conditions 2:2(0). Now if 



fc=i 



ciiyi(t) + ci2X2(0)exp C22(s)ds, t ^ r^. 



(5.5) 



has a unique solution in PC^. 

q 

Lemma 5.2. [18] Suppose a £ PC^: and a = ^ln[ ]^ TT^]' then 

I ^('^fe") = (1 + Cfc)2(Tfc), i = Tfc, 

has a solution z G PC^j if and only if 

/ n ^ cxp[- / a(?/)du]M = 

For the existence of positive periodic solution, we have the following theorem: 

Theorem 5.1. Let us assume (3{t) e PC^^, then there exists a sufficiently small constant cq > 0, such that for each 
P{t) e PCu; satisfying b + a + j < (3{t)S*{t) < & + Q + 7 + C0, there exists a solution {S{t), I{t)) e PC^ x PC^ of 
(2.3) and (2.4) satisfying S{t) < S*{t),I{t) > for aU t. 

Proof. Let ui = S{t) - S*{t), = I{t) in (2.3) and (2.4), then we have 



- p(t)S*(t)u2 - (3it)uiU2, 



t^Tk 



u'2 = [l3{t)S*{t) - (6 + 7 + a)]u2 + I3{t)uiu2, fc = 0, 1, 2 ■ 

. W2(Tjf) = U2(Tfc), 



Denote e{t) = P{t)S*{t) - fi, fi ^ P{t)S*(t), then P{t)S*{t) = 6'(t) + ^ and 6'(t) = 0, thus (5.6) becomes 

u'l ~ —bui — (3{t)S*{t)u2 — P{t)uiU2, t i= Tk 

u'2 = [6'(t) - (fe + 7 + Q!)]w2 + ^"2 + P{t)uiU2, A: 0, 1, 2 • ■ • . 

W2(rj^) = U2{Tk). 



Consider the linear homogeneous system 



-bv^- l3{t)S*{t)v2, t^Tk 

[e{t) - (fe + 7 + a)]w2, A: = 0,1,2. 



^l('^fc ) = (1 ^ Pfc)wi(Tfc), t = Tfe 
, M^k) = ^'2(Tfc). 



(5.6) 



(5.7) 



(5.8) 
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Notice that 

= 9{t) - (6 + 7 + a) = -(6 + 7 + a) 7^ 0, 

cTr=-Mlln[n t4:]>0' 

k=i 

and so (5.8) satisfies the hypotheses of part (1) of Lemma 5.1. Consequently we have the compact linear operator 

given by (5.4). Using L we can equivalently write system (5.7) as the operator 

equation 

(mi, U2) = juL*(wi, U2) + G(mi, U2), (5.9) 

where 

i*(ui,U2) = iLii-mS*it)L2U2),L2U2), 

G{ui,u2) = {Li{-f3{t)S*{t)L2g2iui,u2) + 51(1*1, M2)), £232), 

and (7i(iii,U2) = ~f3(t)uiU2, 52(^1,^2) = f3{t)uiU2. Here L : PCc^ x PC^^ — > ^'C'l^ x -PC'^^ is linear and relatively 
compact, and G : PC^j x PC^j PC^ x PC^ is quasiequicontinuous and relatively compact and satisfies G = o(|j 
(ui, U2) \\pcS) near (0, 0). Note that we have already got the operator equation (5.9) in such a form so that now we can 
apply standard bifurcation theorems and techniques directly to this equation. A non-trivial solution (mi,U2) 7^ (0,0) 
of (5.9) in PGui x PC^j, for some ix E R [R\s the set of the reals), yields a solution (5, /) = (mi + S"*, U2) of the system 
(2.3) and (2.4) for ^ = P{t)S*{t). Solution {S,I) ^ (S'*,0) will be called non-trivial solutions of (2.3) and (2.4). 

To prove Theorem 5.1, we are going to apply some well-known but pretty standard local bifurcation techniques to 
the operator equations (5.9) [17]. The bifurcation can occur only at the non-trivial solutions of the linearized problem 

(W1,W2) =Mi*(«l,^'2), (W1,W2) ^ (0,0), MGi?. (5.10) 

Let (wi,W2) £ PCui X PCui be a solution of equation (5.10) for some iJ. E R, then by the very manner in which L* 
is defined, (ui,W2) satisfies the system 



v{=-bvi- mS*{t)v2, t^Tk 

v'2 = W) + M - (fc + 7 + a)]v2, /c = 0, 1, 2 • • • . 

W2(Tjf) = V2{Tk). 

and conversely. By using Lemma 5.1, we see that (5.11) and hence (5.10) has a non-trivial solution in PC^ x PG^ if 
and only \i 11 = fi* where 

^Ji* = 6 + 7 + a, 

If /i = /i*, then, by part (2) of Lemma 5.1, (5.11) has one independent solution in PG^ x PC^j- Therefore a 
well-known result [17, 20-21] says that bifurcation occurs at this simple eigenvalue; hence there exists a continuum 
[/ = {(//, ui, ^2)} C Rx PG^ X PC^j of non-trivial solutions of (5.9) such that the closure U of U contains (/x*, 0, 0). 
This continuum gives rise to a continuum G = {{fj,; S, I)} C Rx PG^ x PG^ of non-trivial solutions of (2.3) and (2.4) 
whose closure G contains the bifurcation point {fi*;S* ,0). 

To observe that particular solutions in C correspond to solutions {S,I) of (2.3) and (2.4) with the properties 
described in Theorem 3.1, we investigate the nature of the continuum G near the bifurcation point (/i*,0,0) by 
expanding /i and (ui,W2) in the Lyapunov-Schmidt series (for small e) [22]: 

fi = fi* + fiie + • • • , = Uiie + Mi2£^ + • • • ,1 = 1,2 

for Uij e PGui- We substitute these series into the system (5.7) and equate the coefficients of s and e^, we find 

' = -buii - l3{t)S*{t)u2i, ty^Tk 
^ u'2i = [9{t) + n* -{b + j + a)]u2i, fc = 0,l,2.... 

Ull{Tk) = - Pk)uil{Tk), t = Tk 

^ U2liT^) = U2l(Tfe). 
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and 



' "12 = -bui2 - l3{t)S*{t)u22 - /3{t)uiiU21, t ^ Tfc 

"22 = [^(^) + M* - + 7 + a)]u22 + M1U21 + f3{t)uiiU2i, 

Ul2{Tk) = - Pk)ui2iTk), t^Tk 
^ U22iTk) = U22(Tfe), fc 0, 1,2 • ■ • . 



(5.13) 



Obviously (un, U21) G -PC^j x PCui must be a solution of (5.11). We choose the specific solution satisfying the initial 
conditions U2i(0) = 1. Then 



(5.14) 



M2i(t) = exp( / [0{s) + M* - (fe + 7 + a)]ds) > 0. 

JO 

Moreover, unit) is the t^j-periodic solution of the linear equation 



-buii ~ l3{t)S*{t)u2i, ty^Tk 

Ullir^) = (1 -pk)uil{Tk), t Tk 



'■11 



(5.15) 



Hence 



Mil 



Git,s)l3is)S*{s)u2iis)ds, 



Where G{t, s) is Green's function given by 



X{t){l~ X{uj))-yX{s), 0<s<t<uj, 
G{t — ku!, s — joj), juj < s < juj, koj < t < kijj + cj, 



and 



x{t)= n (i-pfeK 

0<rfc<t 



bt 



Note that ui2(t) > for all i, and also since 



X(u:)^\{{l-Pk)e~'- <l, 



k=l 



therefore Green's function G{t,s) > and thus Uii{t) < for all t. Using Lemma 5.2 on the second equation of 
(5.13), we obtain 



(5.16) 



"2i(t)(Mi + f3{t)uii) exp( / [e{s) + /i* - (& + 7 + a)]ds) = 0. 
From the equation of (5.16), we obtain 



^ll = -/3(t)uii(i) > 0. 

Thus we conclude (see [17]) that near the bifurcation point (/x*;0,0) (say, for a sufficiently small constant cq > 0, 
fi* < < fi* + Co ), there is a branch C/+ = {{fi;ui,U2) e U : pi* < ^ < fi* +co,ui{t) < 0,U2{t) > 0}. 

Therefore, we have the piecewise continuous branch C+ of solutions of the from (/i;5, /) G i? x PG^; x PGu), 
satisfying S{t) < S*{t), I{t) > for aU t. Since /x* < < yu* + cq is equivalent to & + a + 7 < (3{t)S*{t) < b + a + j + co, 
which again equivalent to Rq 1+. Hence the proof. 



6. Numerical Simulation 



In our proposed pulse vaccination SIR model with periodic infection rate /3(t) we obtain that the periodic infection- 
free solution is globally stable provided the basic reproductive number Rq is less than 1. On the other hand, periodic 
infection-free solution is unstable if the basic reproductive number is greater than 1 and in this case the disease will 
persist in the population. Moreover, the positive periodic solution exists when Rq —f . 
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In this section we have performed numerical simulation to show the geometric impression of our results. In all 
simulations, the units are set to unity (scaled to unity). 




The time t The time t The proportion ol Ihe susceptible 

(d) (e) (f) 

FIG. 1: Solution of system (2.3) and (2.4) with Ro < 1. Subfigures 1(a), (d), and (b), (e) demonstrate respectively the 
susceptible and infective population with respect to time, and Figure 1 (c), (f) show their corresponding phase-portrait i.e in 
S — I plane. 

To demonstrate the global stability of the system (2.3) and (2.4) we take following set parameter values: b — 0.05, 
K = 1, = 0.25, a = 0.002 and /3(t) = 0.8[1 + 0.2 cos(0.057rt)]. Figure [J depicts the global stability of system (2.3) 
and (2.4) when i?o < 1 with initial conditions S'(O) = 0.408166, /(O) = 0.002, i?(0) = 0.08. Taking a fixed value of 
Tfe = 8 we have varied value the parameter pk- Fig.^a)-(c) shows the solution for pk ~ 0.6, Rq = 0.85909 and for 
Fig.IIId)-(f) we take pk = 0.505, Rq = 0.99892 {Rq 1"). 
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FIG. 2: Unstable infection-free situation of the system (2.3) and (2.4) at Rq > 1. Subfigure 2(a),(d),(g) and (b),(e),(h) 
demonstrate respectively the susceptible and infective population with respect to time t. Subfigure 2(c), (f) and (i) show their 
respective phase portrait i.e in 5 — / plane. 

By choosing Rq > 1 we have shown the unstable infection free solution of the system (2.3) and (2.4) in figure |2] 
Here we have chosen two sets of initial conditions: 5(0) = 0.108166, /(O) = 0.001, R{0) = 0.08 and 5(0) = 0.708166, 
/(O) = 0.2, R{0) = 0.08. The Fig. El;a)-(c) with pk = 0.5 and Rq = 1.00702; The Fig. |3d)-(f) with pk = 0.15 and 
Rq = 1.88779; The Fig. [2Ig)-(i) is the hmiting case with Pk = and Rq = 2.64901. The other parameters have the 
same value as in Fig. 



Now we fix the parameter pk and vary the parameter t^. In figurc|31 the initial values are used 5(0) ~ 0.408166, 
/(O) = 0.002, R{0) = 0.08. Fig. Ei;a)-(c) the parameters are chosen as p = 0.5 and r = 5 with Rq = 0.72972; 
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Fig-EJd)-(f) the parameters are chosen as p = 0.5 and = 13 with Rq = 1.34774, the other parameters as the same 
FigH Fig. O suggests the positive periodic solution is globally asymptotically stable when the i?o > 1- 





FIG. 3: The orbits of the system (2.3) and (2.4)with respect to time series. Subfigure 3(c), (f) depict the orbits on S — / plane, 
and (a), (d) and (b), (e) show the time series for susceptible and infective respectively. 

7. Discussion 



Vaccination strategies are designed and applied to control or eradicate an infection from the population. As for 
many directly transmitted infectious diseases levels of infection oscillate regularly it may be more useful to use a time 
dependent contact rate in place of a constant contact rate. 

The idea of this vaccination strategy is to give the whole susceptible population vaccine at a periodically varying rate 
f3{t): Many vaccination strategies are applied in different parts of the world. Some of these use a periodic vaccination 
strategy. The sort of periodic vaccination rate most commonly used is the pulse periodic function [24] . As the study 
of periodicity and other oscillatory behavior in the incidence of infectious diseases, especially childhood diseases, is 
one of our main goals in this work, we shall study the effect of a general time dependent periodic vaccination strategy 
on the dynamics of these infectious diseases. This new vaccination strategy is designed to vaccinate susceptible of 
all ages using a varying periodic vaccination rate f3{t). This periodic vaccination is intended to keep the disease free 
solution stable by controlling the number of susceptible over the vaccination interval. We have study an SIR epidemic 
model with seasonality in the transmission rate but the results apply equally well to the case where the contact rate 
is constant. 

A lot of infectious diseases have, at least, temporary varying transmission rate due to the seasonal fluctuation or 
changing social behavior. The time dependent epidemic models with pulse vaccination are natural generalization and 
have more applications in real world. The study of the dynamical behaviors of epidemic models with time dependent 
infection rate and pulse vaccination is an important and challenging issue for us. The SIR epidemic model with 
periodic infection rate and pulse vaccination is one of the simple and important epidemic models. In this article, we 
have investigated the dynamical behavior of a classical periodic SIR model with pulse vaccination. We have shown that 
the infection-free periodic solution 0, R*{t)) is globally stable if the basic reproductive number Rq is less than 1, 
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and the infection-free periodic solution {S*{t), 0, R*{t)) is unstable if the basic reproductive number Rq is greater than 
1. In the later case, the disease will uniform persist in the population. Moreover, we have used a standard bifurcation 
technique to show the existence of the positive periodic solutions which arise near the infection-free periodic solution 
{S*{t),0,R*{t)). 

The basic reproductive number is a fundamental parameter and one of the most useful quantities characterize the 
magnitude of an infectious disease transmission [23]. It is the average number of secondary cases that arise from a 
primary case in a fully susceptible population. For the classical SIR model with constant infection rate, the basic 
reproductive number Rq is the product of the average infection period and the average infection in unit time. From 
a deterministic viewpoint, the eradication of an infection requires i?o < 1: Sund the condition for outbreak or for 
maintenance of epidemic infection demands i?o > 1- The basic reproductive number Rq for the SIR model with 
periodic infection rate and pulse vaccination is defined using the average of P{t)S*{t) over one period. It is proved 
that the disease will die out if i?o < Ij and the disease will persist when i?o > 1. The classical SIR model has a 
globally stable endemic equilibrium. The periodic solution of the SIR model with periodic infection rate and pulse 
vaccination is also expected to play the similar role. However, here we have got comparatively few new results on 
the existence of the positive periodic solution of the epidemic models. The bifurcation theorem has been used to 
establish the existence of periodic solution and found that positive periodic solution of our system is globally stable 
when i?o > 1- On the other hand, infection free periodic solution is globally stable if i?o < 1- We have simulated the 
global behaviors of the model for both cases. From the numerical simulation, we conclude that the positive periodic 
solution of the SIR model is unique and globally stable. The uniqueness and stability of the periodic solution is of 
course more challenging problem for the researcher [24] . We expect the uniqueness and global stability of the periodic 
solution will be proved very soon. 
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